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ABSTRACT 

We present a deep color-magnitude diagram for individual stars in the halo of the nearby spiral 
galaxy M81, at a projected distance of 19 kpc, based on data taken with the Advanced Camera for 
Surveys on the Hubble Space Telescope {HST). The CMD reveals a red giant branch that is narrow 
and fairly blue, and a horizontal branch that has stars that lie mostly redward of the RR Lyrae 
instability strip. We derive a mean metallicity of [M/H] = -1.15 ± 0.11 and age of 9 ± 2 Gyr for the 
dominant population in our field, from the shape of the red giant branch, the magnitude of the red 
clump, and the location of the red giant branch bump. We compare our metallicity and age results 
with those found previously for stars in different locations within M81, and in the spheroids of other 
nearby galaxies. 

Subject headings: galaxies :halos, galaxies: stellar content, galaxies: spiral, galaxies:individual(M81) 



1. INTRODUCTION 

Spheroidal stellar populations, such as elliptical galax- 
ies and the halos and bulges of spiral galaxies, contain a 
significant fraction of all stars in the local universe. The 
ages and metallicities of halo stars in these galaxies pro- 
vide some of the most important observational clues to 
the formation and earliest evolution of galaxies, whether 
through studies of globular star clusters (GCs) or, in 
more nearby galaxies, the individual halo stars. Below, 
we broadly summarize what is currently known about the 
ages and metallicities of stars in the spheroids of nearby 
galaxies. 

The spheroid of the Milky Way has two components, a 
centrally concentrated bulge composed of more metal- 
rich stars, and an extended, lower density halo con- 
sisting of more metal-poor stars. The halo of the 
Galaxy shows substructure, with (at least) two chem- 
ically distinct components. Star s of the 'inner halo ' 

have a mean [Fe/Hj; 16 (e.g.. iRvan fc Norrisl [T99TI: 

lAllende Prieto et al.|[200lllvezic et al.ll2008[ ). while stars 
in the 'outer halo' (R>15 k pc) are more meta l-poor, 

with a mean [Fe/H] 2.2 (fCaroUo et al.ll2007[) . The 

halo also contains many stellar streams and tidal debris 
that are probably the remnants of previously accreted 



satellite galaxies (e.g., Ncwberg et al. 2002; Yann v et al 
2003 : Juric et al 200 8; Smith et al. 2009; Ncwber g et al 



20091: iKlement et al.H2609i and references therein). The 



globular clusters in the Galaxy have a weak metallicity 
gradient inside >10 kpc, with no such trend apparent 
beyond 10 kpc out to the most distant cluster located at 
wl20 kpc (Zinn 1985; Harris 2001). The stars and GCs 
in the halo of the Milky Way appe ar to be mostly old , 
with ages between 11 an d 13 Gyr (jCarollo et al.l 120071 : 
iMarin-Franch et al.ll2009tl , indicating that the dominant 



era of accr etion in the Milky W ay took place soon after 
formation (jHammer erall[200l . 

Detailed studies of the stars in M31 suggest that this 
galaxy has had a more active, recent accretion his- 
tory when compared with the Milky Way. M31 has 
a markedl y metal -rich inner spheroid compo nent (e.g., 
T986; Dur rell et all 2001: B rown et alJ 
I2OO6; Gilbert et al.i2007i\that is dom- 



iMould & KristianI 
■2003; Kalirai et al 



inate d by substructu re (Ferguson et al.ll2002tllbata et al.l 
I2OO7I: IRichardson et al.ii200a ). M31 accreted a relatiyely 
mass i ve satellite galaxy o nly ~ 1 Gy r ago ()Font et al.l 
[200I iFardal et alJ [20061 [20071 [20081) . and has had 
other many more minor interactions (jMcConnachie et al.l 
120091 ) as well. Deep images taken with the HST, which 
reach below the main sequence turnoff region for ancient 
stars in M31, indicate the presence of metal-rich, inter- 
mediate age stars (« 6 — 10 Gyr) in fields located at de- 
projected distances of 11, 20, and 35 kpc, with a larger 
contribution from older, more metal-p oor stars increas- 
ing a t larger distances ([Brown et al.l l2003l 120061 . I2007L 

mm- 

Beyond ~ 30 kpc and extending out to at least 
~ 160 kpc, the halo of M31 is do minated by 

(llrwin et al.l 



stars with [Fe/HV 



-1.5 



metal- poor 

20051 iGuh athakurta et"all l2006t iKalirai et al.l 120061 : 



Chapman e t al. 2006), and this component too has copi- 



ous s ubstructure (|Ibata et al.ll2007t iMcConnachie et aD 
|2009| ). Little is known abo ut the ages of these stars 
except for the recent work of iMackev et al.l ([20090 who 
studied the M31 GC MGCl, which is located at a pro- 
jected galactocentric distance of ^^120 kpc. The CMD 
of this cluster reveals a metal abundance of [M/H]~-2.3 
and a horizontal branch morphology that is consistent 
with those of Milky Way GCs at this metal abundance. 
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suggesting a similar age. 

Studies of stars in the spheroids of several more dis- 
tant elliptical galaxies suggest that metal-rich stars dom- 
inate out quite far, probably to at least 10-15-Re//. 
NGC 5128 h a,s a largely metal-r i ch p opulati on at radii 
of 8- 33 kpc ()Harris et all [TQQI [Harris fc HaTrisI [2000l 
|2002|) . Even further out, an analysis of the colors of 
RGB stars, which are sensitive to metallicity, and the 
colors and luminosities of the AGB bump and the RC, 
which are sensitive to both metallicity and age, suggests 
that stars in NGC 5128 at i? ^ 38 kpc have lar gely in- 
term ediate to old ages, with 8 ± 3 Gyr (Rcikub a et alJ 
[2005h . The halos of NGC 3377 and NGC" 3379 have 
a wide range of metallicity, with a significant fraction 
of metal-rich stars, although metal-poor stars begin to 
dominate the spheroid of NGC 3379 beyond w 12Reff. 
The lack of stars brighter than the RGB suggests that 
NGC 3377 and NGC 3379 con tain few young, bright 
AGB stars (IHarris et all I2007allbl) . Like NGC 3379, a 
dominant population of metal-poor stars have been found 
in the spheroid of th e edge-on spiral galaxy NGC 891 
(|Reikuba et alJ[2009l ). 

M81 is an earlier-type spiral (Sab) than either the 
Milky Way or M31. However, studies of the resolved stars 
in the outskirts of M81 have already sh own that it too 
is a co mplex mix of stellar populations. iTikhonov et al.l 
(|2005D used 9 archival WFPC2 images of M81 to show 
that the number density of RGB and AGB stars de- 
creases with galactocentric radius, eventually flattening 
out at a (deprojected distance) of ~ 25 kpc. They at- 
tribute the stars interior to this as part of the thick 
disk, and stars beyond this as part of the halo. CMDs 
of fields interior to 25 kpc reveal relatively young stars 
(jWilliams et al.l 120091 : lOalcanton et al.ll2009t ). and older 
RGB stars with a relatively large range in metallicity 
(-1 < [M/H] < 0; Mouhcine et al. 2005, Wilhams et al. 
2009). These stars are likely locate d in the disk, either 
thin or thick, of M81. More recentlv. [Barker et al.l (|2OO90 
found a flattending of the stellar density profile beyond 
^ 20 kpc (deprojected) from M81, with more metal-poor 
stars ([M/H]^ than in the interior fields, again in- 

dicating that the halo dominates at these distances. The 
only constraint on the ages of stars in the halo of M81 
comes from the measurement of absorption line strengths 
in globular clusters, which suggest ages that are indistin- 
guishable from the GCs in the Galaxy, within the signif- 
icant uncertainties (« 3 — 4 Gyr) (Schroder ct al. 2002). 

In this work, we study the metallicities and ages of 
stars in a field in M81 that lies at a projected distance 
of 19 kpc, using observations taken with the Advanced 
Camera for Surveys on HST. These data are among the 
deepest available for a portion of M81. The rest of this 
paper is organized as follows: Section 2 presents the data 
and reduction; Section 3 presents the CMD for our field 
and discusses general features; Section 4 presents a new 
determination of the distance to M81 from the tip of the 
RGB, and averages this with previous distance determi- 
nations to give our adopted value; Section 5 discusses 
the red clump, and Section 6 presents an estimate of the 
metallicities and ages of stars in our M81 field. Section 7 
presents some new results on a small population of blue 
stars discovered in our field, and Section 8 discusses the 
results in terms of the formation of M81, as well as our 
main conclusions. 



2. OBSERVATIONS AND REDUCTIONS 

Observations of M81 were obtained with the Wide 
Field Channel (WFC) of the Advanced Camera for Sur- 
veys (ACS) on-board the HST, for program GO-10604 
(PI: Sarajedini). A single ACS/WFC pointing cover- 
ing w 3.4' X 3.4' was centered at a = 9''53'"03^20, 
6 = +68°52'03.6" (J2000.0), at a projected distance of 
18' from the center of M81, or 19 kpc at the adopted 
distance of 3.7 Mpc (see Section 4). 

The specific pointing was chosen to lie along the south- 
western semi-minor axis, beyon d the radius sugge sted 
as the edge of M81's thick disk ()Tikhonov et all 12005). 
where the surface brightness of M81 is still high enough 
that we can determine a meaningful metallicity distribu- 
tion from RGB stars. Figure la shows the location of 
our A CS pointing. It also sh ows the WFPC2 field stud- 
ied by iMouchine et al.l (|2005f) , which falls on the edge of 
the suggested thick disk in M81, at a projected distance 
of 12.8' (14 kpc) from the center of M81. Figur e lb shows 
our po inting superposed on HI contours from lYun et al.l 
([1991 : while our field is outside of M81's purported thick 
disk (at a deprojected distance of 30-34 kpc), it still lies 
along the line-of-sight with some of the ma ny HI fila- 
ment s that exist outside M81's optical disk (lYun et all 
fl99l . While some of this HI has been tidally stripped 
from M81 by encounters wit h M82 and NGC 3077, prob- 
ably w 220 - 280 Myr ago (|Yunlll999D . velocities of the 
HI towards our field are consistent with rotation of the 
M81 disk (M.Yun, private communication). 

Our field was observed on 2005 Sept. 11-13 for 5 orbits 
(2 images per orbit) in the F606W filter and 9 orbits (2 
images per orbit) in the F814W filter, using a standard 
4-point dither pattern. Each individual exposure had 
an exposure time of 1247s, leading to a total exposure 
time of 12470s in the F606W fiher, and 22446s in the 
F814W filter. We retrieved the pipeline fiatfielded images 
from the HST archive as well as the pipeline-processed 
multidrizzled combined images from subsets of the data. 

Photometry was performed wit h the ACS mo dule of 
the DOLPHOT software package polphinll2000l ). which 
is specifically designed for point-source photometry of ob- 
jects in the individual FLT images, which were taken di- 
rectly from the HST archive. Object detection and pho- 
tometry was performed on all exposures simultaneously, 
using one of our deep FSIAW drizzled images (derived 
from a subset of 8 images) as the reference frame. The 
DOLPHOT parameters were set to the recommended 
values given in the DOLPHOT User's Guide Q, which 
includes PSF-fitting using pre-derived PSFs, as well as 
the c alculation and application of CTE correction (Riesi 
I2003D and aperture corrections. The instrumental mag- 
nitudes were converted to the VEGAmag HST photo- 
metric system by adopting the revised zeropoints (for 
data taken before July 4, 2006) of 26.420 for F606W and 
25.536 for F814W, which differ by ~ 0.02-0.03 mag from 
the values originally published bv ISirianm "eFaLl (f2005[) . 
We refer to magnitudes in this system by their filter name 
{niFQOGW and mpsiiw), and in cases when we need to 
convert to the Johnson- Cousins systems refer to the fil- 
ters as V and /. 

We determined the photometric completeness limits 

^ http:/ /purcell. as. arizona.edu/dolphot/ 
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and associated errors in our point-source photometry by 
adding and remeasuring 60,000 artificial stars in the ACS 
images. The artificial stars have a range of magnitude 
(24 < mpsiAw < 31) and color (—0.25 < itifqoqw — 
< 2.25). T he Pritchet interpolation function 
(jFleming et al.lfTOQSl ) 



f{m) 



a{m - mum) 



was used to fit the binned completeness fractions / {— 
number of stars detected / number of stars added) as a 
function of magnitude. Here, a is a parameter that mea- 
sures the rate of decline of f(rn) at mum, and mum is 
defined as the 50% completeness limit. To account for 
completeness variations in both color and magnitude, we 
used only the bluest stars {mpGOGW — "mpsiiw = 0) to 
determine the F814W completeness fraction f{F814:W), 
and redder stars {mpeQew — "mpsuw > 1) to de- 
rive / (F606M^). The limiting magnitude (/ = 0.5; 
iHarria 1990) for each filter is mF606W,Hm = 29.73 and 
mpsi4\Y,iim — 28.79. 

3. GENERAL FEATURES IN THE COLOR-MAGNITUDE 
DIAGRAM 

We selected point sources in our field as follows. First, 
we restricted our sample to objects with a DOLPHOT 
category consistent with stars (0BTYPE=1). Of the 
sources that met this criterion, only those that were mea- 
sured on at least 8 of the 10 individual F606W images, 
and on at least 15 of the 18 individual F814W images 
were retained, yielding a list of 13,858 sources over the 
entire ACS field. To remove remaining, slightly resolved 
objects, we applied the following cuts: < 2 (from 
the PSF fitting), and a crowding parameter CROWD 
< 0.4 mag, where the latter rejected objects whose mag- 
nitudes are likely to be affected by neighboring sources. 
Finally, we kept only those objects with sharpness val- 
ues within an empirically-derived hyperbolic function : 
\SHARP\ < 0.01 + 0.16eO-46(-P'8i4W-28)_ of these cri- 
teria were chosen based on visual inspection of sources on 
the ACS images. The mpQ0Qw~'m'F&i4W versus mpsuw 
color-magnitude diagram (CMD) for the 9081 remain- 
ing objects (representing our 'cleaned' sample of stellar 
objects in the M81 field) is shown in the left panel of 
Figure 2, with the CMD of all objects rejected by the 
CHI/SHARP/CROWD cuts plotted in the right hand 
panel. The photometric errors derived from the artificial 
star experiments, as well as the 50% photometric com- 
pleteness levels, are also plotted in the left-hand panel. 

Comparison of the two plots shows that the vast major- 
ity of rejected objects do not correspond to the sequences 
visible in the 'clean' CMD, lending credence to our choice 
of classification criteria. We note that some true stellar 
sources will be removed by our somewhat stringent se- 
lection criteria, but since we are most interested in those 
objects with the best quality photometry, these are ac- 
ceptable 'losses'. 

The cleaned CMD extends from above the tip of the 
first ascent red giant branch (TRGB) to one magnitude 
below the horizontal branch (HB). Stars on the red giant 
branch are older than a few Gyr. The HB has a morphol- 
ogy which is predominantly redward of the RR Lyrae in- 
stability strip, although we cannot unambiguously detect 



an extended blue horizontal branch due to the presence of 
a small population of blue stars {mFeo6W~'n^Fsi4W ~0), 
which will be discussed in Section 7. In fact, the red 
HB resembles a red clump (RC) most commonly seen 
in relatively young (down to « 1 Gyr) and/or metal- 
rich stellar populations. An analysis of the RGB and 
RC stars is presented in Section 5. There is a small but 
noticeable population of stars above the TRGB, likely 
asymptotic giant branch (AGB) stars of intermediate age 
(~5 to 10 Gyr). From the Besangon Galaxy model of 
iRobin et al1(|200l . we estimate that there should be « 7 
(foreground) stars from the Milky Way with 21 < / < 24 
and 1.5 < {V — I) < A in our field of view, so the 
vast majority of stars in our M81 field in this color and 
magnitude range belong to M81. We will explore these 
populations in more detail in the upcoming sections. 

4. DISTANCE TO M81 
4.1. The RGB Luminosity Function 

We determine the luminosity function of the RGB stars 
in our field by defining an 'RGB zone' and thereby ex- 
cluding stars far from the RGB s equence. We use the 
Z^O.OOOl and Z=0.004 models of iGirardi et"all (pOOl 
to provide a general guide to the shape of the RGB in 
both color and luminosity, after first making the models 
redder and bluer by 0.05 mag to accomodate photomet- 
ric uncertainties; see Figure 3. However, we explicitly 
include stars brighter than the tip of the RGB, likely to 
be AGB stars, by allowing a larger range of magnitudes 
than the model predictions. 

In Figure 4, we show the mp^im/ luminosity function 
of the RGB in our field for those objects in the extraction 
region with 22 < 771^^81414' < 28, with a bin width of 
0.08, and uncertainties given by ^/N. We do not make 
a correction for foreground stars, because we found in 
§3 that these are very sparse in this direction. The red 
clump is clearly present at rnF8i4H' 27.8; we will return 
to this feature in Section 5. 

4.2. Tip of the Red Giant Branch 

In order to interpret our CMD further, we need to 
know the distance to M81. Here, we determine the dis- 
tance from the tip of the RGB (TRGB) in the /-band, 
which has proven to be an excellent distan ce indicator for 
metal-poor, old stellar populations (e.g., Lee ct al. 19931 
Sakai ct al. 1996; Ferrarcs c ct al . 2000; BcUazzin i et ^ 
200 It iMcCbnnachie et all [20041: iMouchine et al.l \200^ . 
Visual inspection of our CMD shows a transition in the 
number counts at mpsuw ~ 24.0, albeit with a small 
number of stars. While a small number of stars in the 
upper part of the RGB (there are ^ 110 in our CMD) can 
lead to a possible systematic bi as in the determination 
of It RGB , .Madore fc Freedman (|1995D have shown that 
with at least ~ 100 stars in the top magnitude of the 
RGB (as is the case for our luminosity function), such 
biases are expected to be small. 

To derive a quantitative estimate of the location of 
the RGB tip, we h ave employed the Sobel edge-detection 
filter (S akai et al.l [T996) in order to find the strongest 
discontinuity in the bright part of our LF in Figure 4 
(for the relatively small number of stars in the upper 
RGB, this simple edge-detection method should suffice 
for our purposes). Each object in the LF is modelled as 
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a Gaussian distribution with a dispersion equal to the 
expected photometric error. We use the dispersion in 
the mean of the individual DOLPHOT magnitudes as 
our photometric error. The resulting luminosity function 
$(m) is plotted at the top of Figure 5. To determine the 
location of the RGB tip, we applied to this LF the edge- 
detection algorithm 

i?(m) = $(m + cr(m)) - $(m - cr(m)) (2) 

where a{m) is the mean photo metric error define d from 
our artificial star experiments (jSa kai et al.lll996l ). The 
resultant response function is plotted in the bottom of 
Figure 5, and shows a peak at mj?8i4W — 24.00 ± 0.02, 
where the small formal uncertainty reflects only the 
width of the Sobel function peak. Conve rting to the I- 
band via the synthetic transformations of iSirianni et al.) 
(pOOSi ) (and using the revised zeropoints, as noted above), 
we find Itrgb = 23.97 ± 0.04, where we have adopted 
an uncertainty of 0.03 mag in the transformation from 
to /. The true uncertainty in Itrgb is almost 
certainly larger than that given above, since they do not 
include any uncertainties in the photometry itself. 

The distance to M81 is determined by comparing the 
measured magnitude for the tip of the RGB found above 
with the predicted intrinsic magni tude , Mj , t rg b = 
-4.05 ±0.02 (for [Fe/H]= -1.6) from lRizzi et~all (pOOl . 
This predicted value is very similar t o M i trgb = 
—4.06 ± 0.07 (random error) derived by iFerrarese et al.l 
([2000), and to M,.trgb = -4.04 ± 0.12 (for [Fe/H] = 
-1.7) from IBellazzini etaLl (pOOl . A comparison be- 
tween observed and predicted Mj^trgb magnitudes gives 
an observed distance modulus of {m — M)i — 28.02±0.05 
for M81, or an intrinsic value of {m-M)o = 27.86 ±0.06 
{d — 3.73 ± 0.16 Mpc), where we have assumed a fore- 
ground extinction of Aj = 0.16±0.04 (using E{B — V) = 
0.08;|Schlcgcl ct al. 1998) and adopted a 0.02 uncertainty 
in the foreground reddening to M81. This is consistent 
with determinations from independent TRGB estimates 
and other methods, as discussed below. 

4.3. Previous and Adopted Distance Determinations 

Several different techniques have been used previously 
to estimate the distance to M81, leading to distance mod- 
ulus estimates that differ by approximately 0.1 mag. 

The period-luminosity relationship for Cepheid vari- 
ables is a robust way to de t ermin e the distance to nearby 
galaxies. IFerrarese et al.l (|2OOO0 determined a value of 
(m - M)o = 27.80 ± 0.08 for M81, ba s ed on 30 Cepheids 
discovered in M81 bv lFreedman et al.l (I1994D. usin g ffS'T 
observations. More recently. iMcCommas eFall ([20091 ) 
determined a distance modulus of 27.78 ± 0.05 ± 0.14 
based on 11 Cepheids, in good agreement with the dis- 
tance determined by Ferrarese et al. Here, we adopt the 
McCommas et al. value as the Cepheid based distance to 
M81. 

Estimates of the distance to M81 based on the tip of the 
RGB method tend to give a range of v alues. We found 
(rn-M)o = 27.86 ± 0.06 above, and iTikhonov et~all 
(|2005[) derived a mean value of (m - M)o = 21 m ± 0.04, 
base d on fields which ar e interior to that used here. How- 
ever, iRizzi et all (|2007D derive a much shorter distance 
modulus of (m — M)o = 27.68 ±0.04 from their detection 
of the TRGI3 in an M81 disk field, and .Dalcanton et al.l 



(|2009D derived a value of (m — M)o = 27.77 from deep 
photometry of an ACS field in M81's disk. 

For this work, we adopt a distance of ((m — M)o) = 
27.80 ± 0.08 {d = 3.63 ± 0.14 Mpc) to M81, which is the 
mean value of the five distance determinations described 
above (4 TRGB values and the Cepheid value). The 
quoted error is the dispersion of the input values (each 
of the individual distance determinations have similar 
uncertainties). 

5. LOCATION OF THE RED CLUMP 

Inspection of Figure 2 shows a well-populated horizon- 
tal branch, indicating a population of relatively old stars. 
While our data do not allow for a clean detection of a 
blue horizontal branch, the magnitude and color of the 
red clump (RCQ) allows us to estimate the age and nietal- 
licity of the dominant stellar population in our field, as 
shown in the next section. 

To derive the photometric properties of the RC, we use 
only those stars with colors in the range 0.80 <V — I < 
1.10; a close-up of this part of the CMD is shown in Fig- 
ure 6, along with the luminosity function of stars in this 
color range, where the RC is clearly visible. To derive 
the mean magnitude of this feature, we fit the 'back- 
ground' LF immediately brighter and fainter than the 
RC as a linear function, and fit a single Gaussian to the 
LF once the linear background has been subtracted off. 
We note that the choice of bin size and the exact back- 
ground fit have negligible effect on the mean magnitude 
determined for the RC clump. The resulting fit (and 
la uncertainties) yields Irc = 27.75 ± 0.04. The Gaus- 
sian width of a I = 0.14 ± 0.04, while consistent with 
the expected photometric errors from the artificial star 
experiments, does allow for some intrinsic spread in age 
and/or metallicity. Given our adopted distance modulus 
of ((m - Af)o)= 27.80 ± 0.08, the Lband absolute mag- 
nitude in this region of M81 is Mi{RC) = -0.21 ± 0.10. 
The location of the RC in our dataset is also c onsistent 
with the value derived by (jWilliams et al.l[2009( ) for their 
M81 outer disk fielcfl. 

We note the presence of a small positive feature at / ^ 
27.42 ±0.05 (adopted error based on the bin size used) in 
the LF in Figure 6. This feature is very likely the faint 
signal from the RGB bumjQ, the location on the RGB 

^ Here we equate the RC with the 'red horizontal branch' (RHB), 
where the vast majority of helium-burning stars he redward of the 
instability strip 

^ This may seem odd considering t he large differences in stellar 
population between the two studies. IWilliams et al.l (j^OOSP) quote 
mpgi4\Y,RC = 27.792 ± 0.002 while we get mF8l4W,-BC = 27.78 it 
0.04. The Williams et al. field exhibits a mean age of 2-3 Gyr and a 
metallicity of [M/H]~ —0.7. In our field (see next section) we find 
a mean age of 9 Gyr and a [M/H]= -1.15 ±0.03. Our M81 field is 
~ 6 Gyr older and ^ 0.5 dex more metal-poor than the Williams 
et al. field. A change of —0.5 dex makes the red clump 0.15 mag 
brighter and a change of +6 Gyr makes the red clump 0.17 mag 
fainter (based on Fig. 8 of iWilliams et al.n2009j) . so they offset 
each other resulting in a RC absolute magnitude that is essentially 
unchanged. 

This feature is not the AGB bump, which would be more 
luminous than t he RC by 1 mag (e.g., IReikuba et~all l2005i : 
IWilliams ct al. 2003). We do not see strong evidence for an AGB 
bump here beca use of the relatively low number of stars in our 
field. IWilliams et al.l II2009I ) do see an AGB bump, because their 
field contains many more stars, but do not observe an RGB bump, 
likely because their significantly larger range of metallicities and 
ages washes out this feature. 
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where the hydrogen-burning shell crosses the discontinu- 
ity in chemical ab undances remaining from the stars' con- 
vectiv e envelope (jFerraro et al.lll999t lAlves fc Saraiedinl 
119991 ). The luminosity of the RGB bump depends on 
both the metallicity and the age of the stars, and will 
be used as a consistency check on our results in the next 
section. 

6. MEAN METALLICITY AND AGE OF THE M81 FIELD 

The absolute magnitude of the red clump, as well as the 
shape (slope) and location (color) of the RGB are both 
sensitive to the age and metallicity of the stellar popula- 
tion. However, the RGB is more sensitive to metallicity 
than to age, and the magnitude of the RC is more sensi- 
tive to age than to metallicity, so it is possible to use both 
features together to constrain the age and metallicity of 
the dominant population. 

We comp are the observations with the theoretical 
isochrones of iGirardi et al.l (|2000[ ). assuming scaled-solar 
abundance We first correct the colors and magnitudes 
using the values for reddening and distance noted above. 
Next, we assume an age and use the colors of RGB 
stars with M/<-2.0 to det ermine metalhcity b ased on 
an interpolation within the IGirardi et al.l ()20QClD grid of 
theoretical isochrones with -2.3<[M/H]<0.0. Figure 7 
shows these theoretical RGBs overplotted on our M81 
CMD for an age of 9 Gyr. The Mi <-2.0 magnitude 
range is chosen to minimize the effect of asymptotic 
giant branch stars on the blue side of the RGB (e.g., 
ISaraiedini fc Jablonkal[2005l ). The resultant mean metal- 
licity and the I-band absolute magnitude of the RC de- 
rived above, Mi (RC) = -0.21 zfc 0.10, are used in con- 
junction with the lGirardi et al.l (|200Ql) models for the me- 
dian magnitude of the red clump stars to infer the age of 
the population. This age estimate r is then used to re- 
determine the mean metallicity from the RGB stars and 
the process is repeated. This iterative procedure con- 
verges quickly and yields a mean metallicity of [M/H] 
= —1.15 ±0.03 (random), ±0.11 (systematic), where the 
latter value includes the influence of errors in the distance 
modulus and reddening. The mean age of the dominant 
population is then (r) = 9 ± 2 Gyr, where the quoted un- 
certainty reflects the error in the red clump magnitude, 
which is the dominant source of uncertainty. The resul- 
tant metallicity distribution function (MDF) is shown in 
Fig. 8 wherein the fitted Gaussian profile gives a l-cr 
width of 0.39 dex in [M/H]. 

The dominant signal in the MDF is due to stars around 
Mj ~ —2 to ~ —3. The artificial star tests we have per- 
formed suggest that the mean color error at Mj ~ — 2.4 
is 0.032 mag, which corresponds to <T[m/h] ^ 0.12 dex 
based on the lGirardi et al.l pOOO.) isochrones in the range 
-1.3<[M/H]<-0.7. This suggests a significant fraction of 

^ The determination of stellar metallicities for the RGB stars in 
our field is performed in the HST VEGAmag m,F606W — "^F814VK 
CMD as this makes use of the well-determined photometric zero- 
points for the ACS/WFC instrument. However, the theoretical 
dependence of the red clump magnitude on metallicity and age 
is performed in the ground-based system as it requires knowledge 
of the median I-band RC magnitude for a population of stars, and 
these have only been constructed for the ground-based filters. This 
is unlikely to introduce a significant systematic eff'ect because the 
ACS/WFC F814W passband closely matches its gr ound-based I- 
band counterpart. We also note that the R GBs of thelGirardi et al] 
(|200(]| ) isochrones are identical to those of IMarigo et al.l l|2008I V 



the (T[Af/H] exhibited by the MDF is due to an intrinsic 
range in metallicity among the M81 stars in our field. In 
particular, subtracting the dispersion due to errors from 
the measured dispersion in quadrature yields an intrinsic 
abundance spread of (T[7\////]= 0.37 dex. 

It is also important to keep in mind that we have used 
scaled-solar isochrones for the determination of the MDF 
shown in Fig. 8. The degree to which field stars in 
the M81 halo are enhanced in the a-capture elements 
is currently not known. We can estimate the effect of 
su ch an enhancement on our MDF using the formalism 
of iFerraro et"all (pOOOl ): 

[M/H] = [Fe/H] + logiOmSfa + 0.362) 

where fa represents the factor by which the a-elements 
are enhanced. Thus, if [a/Fe]=+0.3, then /q=2 and 
our computed scaled-solar [M/H] values are larger than 
[Fe/H] by - 0.2 dex. 

We conclude this section by providing a series of con- 
sistency checks on our principle result, which is that the 
mean metallicity and age of stars in our M81 field are 
[M/H]= -1.15 ± 0.11 and 9 ± 2 Gyr, respectively. First, 
the location of the weak RGB bump (described in the 
previous section) at J ~ 27.42 ± 0.05 can be used to pro- 
vide an independent estimate of the age/metallicity of 
the M81 field stars. Indeed, the RGB bump lies above 
the HB/ RC for primarily o l d, more metal-poor pop - 
ulations (jFerraro et all 119991: lATves fc Saraie"dhlil ll99^. 
From our adopted distance modulus, an assumed mean 
color {V - I) = 0.95 ± 0.05 and Ay = 0.27, we find 
MvRaBbumv ^ - 30 ± 0.10. From the relations in 
lAlves fc Saraiedinil ()1999() (and adopting an additional 
0.1 dex uncertainty in those relations), this translates 

to a metalhcity [Fe/H] 1.6 ± 0.2 for an old, 12 Gyr 

population, or a metallicity [Fe/H]^ —1.3 ± 0.2 for a 
9 Gyr population. While the uncertainties in these re- 
sults are large, and the a enhancement unknown, the age 
and metallicity derived from the RGB bump are broadly 
consistent with what we find from the location of the 
RC. A further comparison of our d erived value of Mf 
with predicted values presented bv lReikuba et al.l (I2005D 
their Fig. 22; based on models by iPietrinferni et aLl 



2001 ) are also consistent with a metal-poor population 
with an age of ~ 8 Gyr. 

As an additional check, we compare our CMD with 
that for the Galactic globular cluster NGC 362, stud- 
ied as part of the HST Treasury project GO -10775 
(jSaraiedini et al.l 120071: [Marin-Franch et al.ll2009[ ). This 
cluster has a meta l hcity of [Fe/H]= —1.09 (on the 
[Carretta & Gratton' (1997^ scale), and an age of ~10.5 
Gyr (Marin-Franch et al. 2009), and therefore provides a 
good reference for our field. Figure 9 shows our M81 pho- 
tometry compared with the NGC 362 fiducial sequence 
derived from the HST Treasury project data. The latter 
has been shifted by A(mF8i4w)=13.12 and A{mF606W — 
'7iF8i4w)=0.05. These are based on adopted values of 
(m - M)„ = 14.77 a nd E{B - V) = 0.037 for NGC 
362 (ISaraiedini I I2009t ) and (m - M) Fsiiw = 27.95 and 
E{F606W - F8UW) = 0.080 for M81. We see good 
agreement between the location and shape of the NGC 
362 fiducial and the locus of M81 stars along the red 
clump and the RGB. The color distribution of M81 RGB 
stars is slightly bluer than the NGC 362 RGB fiducial. 
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and is most evident within one magnitude of the M81 
RGB tip. This could be due to NGC 362 being slightly 
more metal-rich (~0.05 dex) and somewhat older (1 to 2 
Gyr) than the mean properties of the M81 field. Taken 
together, these could account for a color difference of 
~0.03 mag in uifgogw — ™_F814VK on the RGB. 

7. BLUE STARS 

Our CMD of the M81 halo field shows a (small) pop- 
ulation of blue stars with mpeoew — ™_F814H' ~ 0, prop- 
erties expected for young 10* yr) stars, but not for 
the stars in the halo of a spiral galaxy. We have checked 
that these are genuine sources with blue colors. They 
are clearly point-like and not obvious background galax- 
ies, and are scattered throughout our ACS field. We 
conclude that these are likely individual, young stars as- 
sociated with MSI. Furthermore, the rather sudden ap- 
pearance of these objects at mpsuw ~ 26 blueward of 
the dominant RGB suggests that these are not RGB or 
RC stars in M81 scattered due to large photometric er- 
rors. An overlay of metal-rich fZ=0 008; [M/H] 0.4) 

isochrones from iMarigo et al.l (|2008f) onto our CMD is 
shown in Fig. 10, and suggests that these blue stars 
have ages between w 200-400 Myr. 

What is the origin of these stars? They have colors 
and magnitudes similar to a stellar population seen in 
the CMD of an out er disk field in M81 (Figure 3 in 
iWilliams et al.|[2009() . Potentially then, there is a sparse 
outer disk in M81, which contributes a very small frac- 
tion of the stars observed in our field. If these stars are 
an extension of the true outer disk of M81, then the disk 
plane extends out ^34-40 kpc from M81. 

A second possibility is that these stars formed in the 
disk of M81 or in one if its companions, and then were 
carried to this location along with the HI tidal debris 
that is coincident with our field, likely as a result of the 
most recent galactic interactions between M81, M82 and 
NGC 3077, which occurred ~ 220 - 280 Myr ago (|Ym] 
119991 see also Weisz et al. 2008). 

The observed HI gas along the line-of-sight to our ACS 
field has a colum n density of only Nhi ~ 10^°cm^^ 
(|Yun et al.|[T99l . Recent (~ 30 - 70 Myr) star for- 
mation has been found not only in the outskirts of 
M81 itself (|Barker et al.ll2009f ) but al so in nearby tidal 
dwarf galaxies su ch as Ho lmberg IX dWeisz et al.l 120081 : 
iSabbi et al.|[2008l) . BK3N (|Makarova et al."2002^, and in 
other nearby tidal debris yP urrcll et al. 2004; Davidgg 
[20081 IChiboucas et~all [2001 iMouhcine fc Ibatal 1201(1 7 
However, in all of these cases, the recent star forma- 
tion has occurred in areas with higher HI column den- 
sities {Nhi > 4 — 8 X 10^°cm~^) than that observed 
in our field, and is thus consistent with the lack of re- 
cent (i.e. T < 10* yr) star formation i n our field. This 
is also consistent with the findings of iMavbhate et al.l 
()2007|) . who find recent low- level star formation at lev- 
els Nhi ^ 4 X lO^^cm"^). We note that while recent 
episodes of star formation in the outermost regions of 
late types galaxies have been observed at column den- 
sities similar to that in our field (^ 10^ °cm~^; e.g. 
Ide Blok fc WaIteI1[2003t iThilker et al1[200l . we do not 
see evidence of this in our (small) field. 

We conclude that while we do see evidence for blue 
stars in our field, these make only a small contribution 
to the stellar population in our M81 field. 



8. DISCUSSION AND CONCLUSIONS 

The main result from our analysis of a field at a (pro- 
jected) distance of 18' (19 kpc) from the center of M81, 
just beyond the suggested edge of the thick disk, is that 
the dominant population has a mean age of 9±2 Gyr and 
a peak metallicity of [M/H]= — 1.15± 0.11. Below, we 
discuss these results in the context of the formation his- 
tory of M81, and compare with results from the spheroids 
of other nearby galaxies. 

The color distribution of stars on the RGB in our M81 
field is relatively blue and rather narrow. We derived the 
metallicity distribution of these RGB stars by comparing 
with theoretical isochrones, and find a mean [M/H] of - 
1.15, with an intrinsic dispersion of 0.37 dex as shown 
in Figure 8. Our field is dominated by stars that are 
more metal-poor than stars located closer to the cen- 
ter of M81; indeed, our ACS field contained the most 
metal-poor stars in any M81 field observed thus far. In 
Figure 11, we compare our MDF (for an assumed age o f 
9 Gyr) with that determined by iMouchine et al.l ()2005l ) 
(for globular cluster-like ages; ~ 12 Gyr for a field that 
is also along the southwest minor axis, but somewhat 
closer to the center of M81 with a projected distance 
of ^14 kpc (see Figure 1). While the general shape for 
the MDF is similar in both fields, the Mouhcine et al. 
field is clearly dominated by stars which are more metal- 
rich by sa 0.4 dex than the stars studied here (where 
we have already accounted for the differences in the as- 
sumed ages of the MDFs - adopting a slightly older (by 

3 Gyr) age would effectively decrease the mean metal- 
licity of our MDF by ~ 0.1 dex; making the difference 
between the MDFs larger). This metallicity difference is 
not likely to be due to uncertainties in the colors or in 
the adopted reddening values; to match a 0.4 dex differ- 
ence in [M/H], our mean RGB would need to be redder 
by ^ 0.1 in (mp a Qaw — t tt-fhuw)- Simila r ly, the MDFs of 
IMouchine etall (j20Q5D ; IMouhcine et al.l ([2006') are more 
m etal-poor than the out er (thick?) disk field presented 
bv IWilliams et al.l (|2009l ). a field dominated by relatively 
old (2/3 of stars formed more than ^ 8 Gyr ago) but 
relatively metal-rich stars. 

This sequence of declining metallicity with location in- 
dicates that either there is a strong metallicity gradient in 
the disk of M81, or that our field (as well as the Mouhcine 
et al. field ) samples a diff e rent s tructural component of 
M81 than IWilliams et al.l (j2009f ). The markedly lower 
metallicity between our field and that of the outer disk 
of M81, combined with the relatively old ages inferred for 
these stars (discussed below) and the lack of (much) re- 
cent star formation, provide tentative evidence that our 
field is largely dominated by stars in the halo of M81, 
rather than in a thick disk. 

In a recent wide-field study of bright stars in M81's 
outer regions [Barker et al.l (|2009D suggest the presence 
of a faint, structurally distinct component surrounding 
M81 with a mean metalhcity [M/H]~ -1.1 ± 0.3; the 
large uncertainty is due to the rather limited part of the 
RGB observable in their ground-based data. This metal- 
licity is slightly higher than (but consistent with) the 
value found here [3. Our field is not located within the 

^ it is important to note that [Barker et aLI II2009I ) assumed an 
age of 10 Gyr; if instead they had assumed 9 Gyr as we derive here, 
this would increase their derived metallicity less than 0.1 dex 
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region surveyed bv lBarker et al.l (|2009[ ). but lies at a de- 
projected distance of 34-40 kpc from M81, comparable 
to the distance of the fields used for their metallicity es- 
timate. This extended component has properties similar 
(but not identical) to either the Milky Way thick disk 
or the halo, making the exact nature of this structure 
somewhat uncertain. Future deep observations of M81's 
outer regions will be needed to help establish the true 
nature of this feature. 

We can compare the MDF that we derive for our 
20 kpc field in M81 with the metallicities found for 
the spheroids/halos of other nearby galaxies. At galac- 
tocentric radii of i? < 20 — 30 kpc, similar to the 
location of our M81 pointing, the spheroid of M31 
is dominated by a largely metal-rich ( [Fe/H]^ 
population (e.g., [Pur rell et al. 2001; Kal irai et al.l 1200(1 : 
iGilbert et alll2007nRichardson et al.ll2008l ). The stars in 
the halo of the Milky Way are metal-poor, with typ- 
ical [Fe/H] values of —1.4 to —1.6 for the inner halo 
(out to fa 10 - 15 kpc) (|Rvan fc N orris"1991': llvezi( 5 et al.l 
120081) and « -2.2 beyond this (Carollo et al. 200^ and 
references within). The halo of the E/SO galaxy NGC 
3379 (at a galactocentric distance of 33 kpc) shows a 
broad metallicity distribution with a significant numbe r 
of metal-poor ([M/H]< -0.7) stars (|Harris et al.ll2007bD . 
In contrast, the large elliptical galaxy NGC 5128 has 
a metal-rich halo with [M/H]= -0.6 to - 0.7 at all lo- 
cations studied thus far, from 8-38 kpc (iHarris et al.l 
[199^, 'Harns_&Harris "200^, [200l IR.eikuba et al.ll2005n . 
Similarly, the elliptical galaxy NGC 3377 also shows 
a lar gely metal-rich popu lation at small galactocentric 
radii iHarris et al.l l2007aD . 

Taken together, large, nearby galaxies exhibit a wide 
diversity of halo metallicities, showing no clear trends 
with parent galaxy mass or morphological type (although 
there are suggesti ons that halo metallicit y increases with 
galaxy mass, e.g.. iMouhcine et "al]|2006[ ). Indeed, it ap- 
pears that variations in the metallicity of halo stars even 
at different locations within a given 'halo ' could make any 
trends with parent galaxy properties difficult to quantify. 
For example, the metallicity of stars in our single M81 
field are similar to that of the 'halo' of NGC 3379, the 
'outer halo' of M31, and the 'inner halo' of the Milky 
Way (assuming [a/Fe] ^ +0.3 — 0.4). If we have indeed 
sampled the halo of M81, an open question is whether an 
even more metal-poor population exists at larger galac- 
tocentric radii; again, future studies of more distant M81 
field would be worthwhile. 

While our CMD does not reach the main sequence 
turnoff region, the most reliable age-sensitive feature in 
the CMD, it still holds clues to the dominant age of the 
stars in this region of M81. Formal fits to the RC, RGB, 
and RGB bump give consistent estimates for the mean 
age of the stars of 9 ± 2 Gyr. The fact that we do not ob- 
serve a strong population of blue horizontal branch stars 
is also consistent with these results. 

We can compare the ages estimated for stars in the 
spheroids of different galaxies, although it should be 
noted that these estimates come from different tech- 
niquefl The dominant population of our M81 field is 

''' while different studies assume differing ages (e.g. ,12 Gyr or 
14 Gyr) for the 'old' populations, the effect of these differences on 
metallicity are small, and do not affect our general conclusions 



slightly younger but consistent with the mean ages (9.7, 
11.0, and 10.5 Gyr, for fields at distances of 11, 21, 
an d 35 kpc from M31 ) of the M31 halo fields studied 
bv I Brown et al.l ()2008L and references within). The halo 
of M31 also harbors globular clusters as old as th ose in 
the Milky Way (Alves-Brit o et al...2009,: Ma et al.i r2009): 
however, the presence or absence of any radial age trends 
among this population is still largely an open question. 

In contrast, the stars in our M81 field appear to be 
slightly younger (by « 2 — 3 Gyr, but just within the age 
uncer tainties), than the oldest stars in the Milky Way 
halo (IM arm- Franch et al.1l2009h . as our comparison with 
the globular cluster NGC 362 illustrated above. While 
there is an age spread of a few Gyr among the Milky 
Way GCs with the inner halo cluste rs being older in th e 
mean than those of the outer halo (|Dotter et al.ll2069| ). 
the ages of globular clusters in M81 (estimated from 
integrated spectra) are indistinguishable from Galac- 
tic globular clusters (within the sizeable uncertainties 
iSchroder et al.ll2002l ). It should be noted that these clus- 
ters are all projected at smaller distances in M81 than 
the location of our field. We note that our data do not 
rule out the presence of a minority population of older 
stars, as ancient as the majority Galactic globular cluster 
system, in this region of M81. 

In the case of NGC 5128. lReikuba et al.l (|2005l ) investi- 
gated the ages of field stars at a distance of 38 kpc using 
the magnitude of the RC, and found that it is consistent 
with an old population with a mean age of 8 ± 3 Gyr, 
very similar to the ages that we observe in the M81 halo. 

Our age and metallicity results for halo stars in 
M81 are in broad agreem ent with predictions from 
some ACDM m odels (e.g., iBullock fc JohnstonI [20051: 
iFont et al.lf2008f ). where the outer (i? > 10 - 20 kpc) 
halos of L* galaxies are expected to be formed from 
the accreted (over a relatively long timescale) debris of 
somewhat lower-mass progenitor galaxies than the more 
metal-rich systems that (for instance) resulted in the for- 
mation of the more m etal-rich Giant Southern Stream 
in M31 (Fard al et al.l [2"006. 2008; Font ct al. 2008). In 
this scenario, the outer halos of luminous galaxies like 
M81 will have more metal-poor and somewhat younger 
combined stellar populations than that of the (slightly) 
older, more metal-rich inner halo /spheroid. Our results 
are broadly consistent with these predictions. 

9. SUMMARY 

We have presented deep photometry of a field located 
20 kpc from the center of M81, believed to be dominated 
by halo stars. The resulting CMD, based on HST/ACS 
imaging, reveals a relatively blue RGB and a HB popu- 
lated predominantly redward of the instability strip (i.e., 
a red clump). From the shape of the RGB, the mag- 
nitude of the red clump, and the location of the RGB 
bump, we derive a mean metallicity of [M/H] — -1.15 
±0.11 and an age of 9 ± 2 Gyr for the dominant pop- 
ulation in this region of M81's halo. This is the lowest 
metallicity found for stars in any portion of M81, almost 
certainly because our field is more distant than any other 
portion of M81 studied to date. Our age and metallic- 
ity results are broadly consistent with the predictions 
from ACDM models, which predict lower metallicities 
and younger mean ages for outer halos of spiral galax- 
ies than for inner halos. Future studies that probe the 
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metallicities and ages of stars even further out in the 
halo of M81 will allow for better constraints on the early 
formation of this galaxy. 
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Figure 1. (Left:) The location and field of view of our ACS field ( at lower right), s uperimposed on a red Digitized Sky Survey image of 
M81. The ACS field is located somewhat beyond the WFPC2 field of'Mouchine^et_aIl! (2005), which is also shown. The ellipse indicates the 
22' (~ 25 kpc) extent of M 81's thick disk a s suggested by Tikhonov ct al. (2005), inclined at 59°. The box at upper center is the location 
of the deep ACS study by IWilliams et al.l (2009); Dalcanton et al. (2009) (Right:) The same thick disk and ACS field pointings as the 
figure on the left are shown, overplotted on the greyscale HI map of the M81 system bv lYun et al.l HI 9941 ). where the contours represent 
N{HI) = 0.5, 1.0, 2.5, 5.0, 7.5 and 10 x lO'^^cm"^. While our ACS field is superposed on some of the HI in the system, the column densities 
at this location are lower than that expected in order for star formation to occur. Both images are roughly 48' high and 57' wide; north is 
up, and east is to the left. 
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Figure 2. (Left:) The _F814iy, F606W — FSUW color-magnitude diagram for all stellar sources measured in our ACS field, after removal 
of objects which did not meet the selection criteria (x'^, sharpness, and crowding) as defined in the text. The error bars on the left denote 
the average errors in the photometry based on the artificial star experiments. The CMD is dominated by an old red giant branch as well 
as a horizontal branch/red clump at mpfn^w ^ 28. A sequence of blue (mi?go6W ~ ^FSitw ^ 0) objects is also seen extending upward to 
^FSltW ~ 26. (Right:) The CMD of all objects measured by DOLPHOT that were rejected by our selection criteria. The dashed lines 
in each figure denote the 50% completeness limits in both mj^goew and 'tif814VK- 
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F606W - F814W 

Figure 3. M81 halo CMD showing the loci used to define our RGB sample. The 12.6 Gyr Girardi et al. (2002) Z=0.0001 (left) and 
Z=0.004 (right) models have been shifted outward by 0.05 magnitude to account for photometric spread. The models were also shifted via 
a distance modulus (m - M)f8UW = 27.95, and reddened by E(F606W - F8UW) = 0.08. 
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Figure 4. 01^81414'' luminosity function for all stars between the RGB loci in Figure 4. No corrections for photometric incompleteness 
have been made, nor have any (small) foreground star contributions been removed. The uncertainties are simply the v'JV statistics for the 
objects detected in each 0.075 mag bin. The dashed line at mpgnw = 28.79 denotes the 50% completeness level of the photometry. The 
solid line at mpsuw = 24.00 represents the location of the RGB tip. 
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Figure 5. The top panel shows a smoothed luminosity function formed by co-adding the Gaussian representation of stars in our RGB 
subsample (see text). The bottom panel is the result of applying the Sobel edge-detection algorithm. The peaJc at mpsi4WTRGB = 
24.00 it 0.02 (internal error only) is the most likely location of the RGB tip. 



14 



Durrell et al. 




Figure 6. (Left:) A close up of the region around the red clump and horizontal branch in the I, V — I color-magnitude diagram 
(transformed from the ACS magnitudes using the Sirianni et al. (2005) transformations - see text). The box illustrates the region within 
which the red clump magnitudes have been derived. (Right:) The binned I band luminosity function for objects within the boxed region. 
The LF has not been corrected for photometric incompleteness. The solid line shows the combined linear background and best-fit Gaussian 
to the location of the red clump, at / = 27.76 it 0.04, and a dispersion a = 0.12Tq'o4 
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Figure 7. Our M81 field CMD overplotted with tiie 9 Gyr RGB models from Marigo et al.' (2008'), for metallicities of (from left to right) 
[M/H]= -2.3, -1.7, -1.3, -0.7, -0.4, 0.0. All models have been shifted by E{F606W - F814W) = 0.08 and (m - M)f814VK = 27.95. 
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Figure 8. A plot of the derived metallicity distribution function for stars in ou r M81 field. The solid line is the binned histogram 
of metallicity determinations for RGB stars with Af/<— 2.0 using the Girardi et al.l II2000I ) isochrones (see text). The dashed line is the 
Gaussian fit to the core of this distribution. We assume an age of log (r/yr) = 9.95. 
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Figure 9. The CMD of our M81 halo field, overplotted with the fiducial of the Milky Way globular cluster NGC 362 from the HST Treasury 
project, and the bounding box of the red IfB. The fiducial has been shifted by A(mi^606W ~^F8l4w) = 0-05 a^iid A"^F8141V ~ = 13.12 
to match the distance and reddening of M81. 




Figure 10. Our M81 field CMD overplottcd with Z=0.008 isochrones from lMarigo et al.l H2008I ). The ages are (from top to bottom) 100 
Myr, 200 Myr, 400 Myr. As before, all models have been shifted by E(F606VF - F814W) = 0.08 and (m - M)fsi4.w = 27.95. 
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Figur e 11. Comparison between our MDF (for an age of 9 Gyr) [solid line] and that of the MDF of an interior field from IMouchine et al.l 
II2005I) (for an assumed age of ^ 12 Gyr; dashed line) 



